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'^ . Abstract. Study of nucleosynthesis in accretion disks around black holes was ini- 

tiated by Chakrabarti et al. (1987). In the present work we do the similar analysis 
^ ■ using the state-of-the-art disk model, namely, Advective Accretion Disks. Dur- 

^^ , ing the infall, matter temperature and density are generally increased which are 

ji^ ' first computed. These quantities are used to obtain local changes in composition, 

f— ^ ■ amount of nuclear energy released or absorbed, etc. under various inflow condi- 

tions. In the cases where the magnetic viscosity is dominant neutron torus may 
2^ ' be formed. We also talk about the fate of Li"^ and D during the accretion. The 

outflowing winds from the disk could carry the new isotopes produced by nucle- 
osynthesis and contaminate the surroundings. From the degree of contamination, 
one could pinpoint the inflow parameters. 
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^ , 1. Introduction 



There are many observational evidences where the incoming matter has 
the potential to become as hot as its virial temperature T^iriai ~ 10^^ K 
[1]. Through various cooling effects, this incoming matter is usually cooled 
down to produce hard and soft states [2]. In the accretion disk, matter 
in the sub-Keplerian region generally remains hotter than Keplerian disks. 
The matter is so hot that after big-bang nucleosynthesis this is the most 
fevourable temperature to produce significant nuclear reactions. The energy 
generation due to nucleosynthesis could be high enough to destabilize the 
flow and the modified composition may come out through winds to affect 
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the nietallicity of the galaxy [3-7]. Previous works on nucleosynthesis in disk 
was done for cooler thick accretion disks. Since the sub-Keplerian region is 
much hotter than of Keplerian region and also than the central temperature 
(~ lO^K) of stars, presently we are interested to study nucleosynthesis in 
hot sub-Keplerian region of accretion disks. 

2. Basic Equations and Physical Systems 

In 1981 Paczyhski & Bisnovatyi-Kogan [8] initiated the study of viscous 
transonic flow although the global solutions of advective accretion disks 
were obtained much later [9] which we use here. In the advective disks, 
matter must have radial motion which is transonic. The supersonic flow 
must be sub-Keplerian and therefore must deviate from a Keplerian disk 
away from the black hole. The basic equations which matter obeys while 
falling towards the black hole from the boundary between Keplerian and 
sub-Keplerian region are given below (for details, see, [9]): 

(a) The radial momentum equation: 

ax p ax x-^ 

(b) The continuity equation: 

-^(Sxi9) = 0, (16) 

dx 

(c) The azimuthal momentum equation: 

„dX(x) 1 d , 9^^, , ,, , 

(d) The entropy equation: 

2napTDh{x) da a^'dh{x) dp 



J dx 7 dx 



/Q+ (W) 



where the equation of state is chosen as a^ = ^^. Here, A is the specific 
angular momentum of the infalling matter, Xxep is that in the Keplerian 
region is defined as A|^ = 2(x-i)'-^ t-''^]' ^ ^® vertically integrated density, 
Wx(j) is the stress tensor, a is the sound speed and h{x) is the half thickness 
of the disk (~ ax^''^{x — 1)), n = ^— j- is the polytropic index, / is the 
cooling factor which is kept constant throughout our study, Q^ is the heat 
generation due to the viscous effect of the disk. For the time being we are 
neglecting the magnetic heating term. 
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During infall different nuclear reactions take place and nuclear energy 
is released. Here, our study is exploratory so in the heating term Q+, we do 
not include the heating due to nuclear reactions. (Work including nuclear 
energy release term is in [6].) Another parameter /? is defined as ratio of 
gas pressure to total pressure, which is assumed to be a constant value 
throughout a particular case. Actually, the factor /? is used to take into 
account the cooling due to Comptonization. To compute the temperature 
of the Comptonized flow in the advective region which may or may not 
have shocks, we follow Chakrabarti & Titarchuk [2] and Chakrabarti's [11] 
works and method. The temperature is computed from. 

It is seen that due to hotter nature of the advective disk especially when 
accretion rate is low, Compton cooling is negligible, the major precess of 
hydrogen burning is the rapid proton capture process, which operates at 
T ^ 0.5 X lO^K which is much higher than the operating temperature of 
PP chain (operates at T ~ 0.01 - 0.2 x lO'^K) and CNO cycle (operates at 
T ~ 0.02 — 0.5 X lO^K) which take place in the case of stellar nucleosynthesis 
where temperature is much lower. Also in stellar case, in different radii same 
sets of reaction take place but in the case of disk, in different radii different 
reactions (or different sets of reaction) can take place simultaneously. These 
are the basic differences between the nucleosynthesis in stars and disks. 

For simplicity, we take the solar abundance as the initial abundance of 
the disk and our computation starts where matter leaves a Keplerian disk. 
According to [2] and [11], the black hole remains in hard states when vis- 
cosity and accretion rate are smaller. In this case, xk (at radius xk matter 
deviates from Keplerian to sub-Keplerian region) is large. In this parame- 
ter range the protons remain hot {Tp ~ 1 — 10 x lO^K). The corresponding 
factor /(= 1 — Q^ /Q~) is not low enough to cool down the disk, (in [1], 
it is indicated that rajo? is a good indication of the cooling efficiency of 
the hot fiow), where Q^ and Q~ are the heat gain and heat loss due to 
viscosity of the disk. 

We have studied a large region of parameter space with 0.0001 ,^ a ,^ 1, 
0.001 < m < 100, 0.01 < /? < 1, 4/3 < 7 < 5/3. We study a case with a 
standing shock as well. In selecting the reaction network we kept in mind 
the fact that hotter fiows may produce heavier elements through triple- 
a and rapid proton and a capture processes. Furthermore due to photo- 
dissociation significant neutrons may be produced and there is a possibility 
of production of neutron rich isotopes. Thus, we consider sufficient number 
of isotopes on either side of the stability line. The network thus contains 
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protons, neutrons, till ^^Ge - altogether 255 nuclear species. The standard 
reaction rates were taken [6]. 

3. Results 

Here we present a typical case containing a shock wave in the advective 
region [6]. We express the length in the unit of one Schwarzschild radius 
which is -2 where M is the mass of the black hole, velocity is expressed 
in the unit of velocity of light c and the unit of time is 3 . We use the 
mass of the black hole M/Mq = 10 {Mq = solar mass), Il-stress viscosity 
parameter an = 0.05, the location of the inner sonic point Xj„ = 2.8695, 
the value of the specific angular momentum at the inner edge of the black 
hole Xin = 1.6, the polytropic index 7 = 4/3 as free parameters. The net 
accretion rate rh = 1 in the unit of Eddington rate, cooling factor due 
to Comptonization /3 = 0.03, xk = 481. The proton temperature (in the 
unit of 10^), velocity distribution (in the units of 10^*^ cm sec~^), density 
distribution (in the unit of 2 x 10^^ gm cm^^) are shown in Fig. 1(a). 

In Fig. lb, we show composition changes close to the black hole both 
for the shock-free branch (dotted curves) and the shocked branch of the 
solution (solid curves). Only prominent elements are plotted. The difference 
between the shocked and the shock- free cases is that, in the shock case the 
similar burning takes place farther away from the black hole because of 
much higher temperature in the post-shock region. A significant amount of 
the neutron (with a final abundance of 1^ '^ 10~"^) is produced due to photo- 
dissociation process. Note that closer to the black hole, ^\J, ^^O, '^'^Mg and 
'^^Si are all destroyed completely. Among the new species which are formed 
closer to the black hole are ^^Si, ^^i, ^^Cr. Note that the final abundance 
of ^''A^e is significantly higher than the initial value. Thus a significant 
metallicity could be supplied by winds from the centrifugal barrier. In Fig. 
Ic we show the change of abundance of neutron (n), deuterium {D) and 
lithium (jLi). It is noted that near black hole a significant amount of neutron 
is formed although initially neutron abundance was almost zero. Also D and 
"^Li are totally burnt out near black hole which is against the major claim 
of Yi & Narayan [13] which found significant lithium in the disk. It is true 
that due to spallation reaction, i.e., 

^He+^He^'^ Li+p 

'^Li may be formed in the disk but due to photo-dissociation in high tem- 
perature all ^if e are burnt out before forming "^Li i.e. the formation rate of 
'^He from D is much slower than the burning rate of it. Yi & Narayan [13] 
do not include the possibility of photo-dissociation in the hot disk. 



NUCLEOSYNTHESIS AROUND BLACK HOLES 







1 1 1 M 1 1 






5 
4 


r\ 


\ n 


P-a 


(a)^ 


3 




V\\ 




- 




- 


■A \ 




- 


2 


j 




\ 


— 


1 


r\ 


1 1 1 1 1 1 


1 1 1 1 1 1^ 


T-f-r-i 1 i-i 1. 



O 



-3 







_, , , , 1 , , , , 1 , , , , 1 , , , , 1 , 


(d); 


10 


- 






r-~^ 


-_____; 





- 






Q„u= 


- 


10 


- 


/ 


\ / 




- 


20 


/y "^ 


- 




- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 r 



0.4 0.6 

log(x) 



O 0.5 1 1.5 2 2.5 
log(x) 



Figure 1 : Variation of (a) proton temperature (Tg), radial velocity uio and den- 
sity distribution /9_8 (b) matter abundance Y in logarithmic scale (c) neutron, 
deuterium and lithium abundance Y in logarithmic scale and (d) nuclear energy 
release and absorption as a functions of logarithmic radial distance x. See text for 
parameters. Solutions in the stable branch with shocks are solid curves and those 
without the shock are dotted in (a-d). At the shock, temperature and density 
rise and velocity lower significantly and cause a significant change in abundance 
even farther out. Shock induced winds may cause substantial contamination of the 
galactic composition when parameters are chosen from these regions [6] . 



In Fig. Id, we show nuclear energy release/absorption for the flow in 
in units of erg sec^^ gni^^. Solid curve represents the nuclear energy re- 
lease/absorption for the shocked flow and the dotted curve is that for unsta- 
ble shock-free flow. As matter leaves the Keplerian region, the rapid proton 
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Figure 2 : The convergence of the neutron abundance through successive itera- 
tions in a very hot advective disk. From bottom to top curves 1st, 4th, 7th and 
11th iteration results are shown. A neutron torus with a significant abundance is 
formed in this case [15]. 



capture such as, p +^^0 — >^^ A^ +^He etc., burn hydrogen and releases 
energy to the disk. At around x = 50, Z) — s- n -|-p dissociates D and the en- 
dothermic reaction causes the nuclear energy release to become 'negative', 
i.e., a huge amount of energy is absorbed from the disk. At around x = 15 
the energy release is again dominated by the original processes because no 
deuterium is left to burn. Due to excessive temperature, immediately ^He 
breaks down into deuterium and through dissociation of D again a huge 
amount of energy is absorbed from the disk. It is noted that energy ab- 
sorption due to photo-dissociation as well as the magnitude of the energy 
release due to proton capture process and that due to viscous dissipation 
{Q~^) are very similar (save the region where endothermic reactions domi- 
nate). This suggests that even with nuclear reactions, at least some part of 
the advective disk may be perfectly stable. 

We now present another interesting case where lower accretion rate 
(m = 0.01) but higher viscosity (0.2) were used and the efficiency of cooling 
is not 100% (/ = 0.1). That means that the temperature of the flow is high 
(/9 = 0.1, maximum temperature Tg""^^ = 11). In this case xk = 8.8, if the 
high viscosity is due to stochastic magnetic field, protons would be drifted 
towards the black hole due to magnetic viscosity, but the neutrons will 
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not be drifted [13] till they decay. This principle has been used to do the 
simulation in this case. The modified composition in one sweep is allowed 
to interact with freshly accreting matter with the understanding that the 
accumulated neutrons do not drift radially. After few iterations or sweeps 
the steady distribution of the composition may be achieved. Figure 2 shows 
the neutron distributions in iteration numbers 1, 4, 7 &: 11 respectively 
(from bottom to top curves) in the advective region. The formation of a 
'neutron torus' is very apparent in this result and generally in all the hot 
advective flows. In 1987 Hogan &: Applegate [14] showed that formation of 
neutron torus is possible with high accretion rate. But high accretion rate 
means high rate of photon to dump into sub-Keplerian region and high 
rate of inverse Compton process through which matter cool down, that is 
why photo-dissociation will be less prominent. Also formation of neutron is 
possible through the photo-dissociation of deuterium in the hot disk which 
is physically possible prominently in our parameter region, where neutron 
torus is formed. Details are in Chakrabarti & Mukhopadhyay [15]. 

4. Discussions and Conclusions 

In this paper, we have explored the possibility of nuclear reactions in ad- 
vective accretion flows around black holes. Temperature in this region is 
controlled by the efficiencies of bremsstrahlung and Comptonization pro- 
cesses [2, 7]. For a higher Keplerian rate and higher viscosity, the inner edge 
of the Keplerian component comes closer to the black hole and the advec- 
tive region becomes cooler [2, 9]. However, as the viscosity is decreased, 
the inner edge of the Keplerian component moves away and the Compton 
cooling becomes less efficient. 

The composition changes especially in the centrifugal pressure sup- 
ported denser region, where matter is hotter and slowly moving. Since 
centrifugal pressure supported region can be treated as an effective surface 
of the black hole which may generate winds and outflows in the same way 
as the stellar surface, one could envisage that the winds produced in this 
region would carry away modified composition [16-18]. In very hot disks, a 
significant amount of free neutrons are produced which, while coming out 
through winds may recombine with outflowing protons at a cooler envi- 
ronment to possibly form deuteriums. A few related questions have been 
asked lately: Can lithium in the universe be produced in black hole accre- 
tion [12,19]? We believe that this is not possible. When the full network is 
used we find that the hotter disks where spallation would have been im- 
portant also heliums photo-dissociate into deuteriums and then to protons 
and neutrons before any significant production of lithiums. Another ques- 
tion is: Could the metallicity of the galaxy be explained, at least partially. 
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by nuclear reactions? We believe that this is quite possible. Details are in 
[6]. 

Another important thing which we find that in the case of hot inflows 
formation of neutron tori is a very distinct possibility [15]. Presence of a 
neutron torus around a black hole would help the formation of neutron rich 
species as well, a process hitherto attributed to the supernovae explosions 
only. It can also help production of Li on the companion star surface (see 
[6] and references therein). 

The advective disks as we know today do not perfectly match with a 
Keplerian disk. The shear, i.e., dVt/dx is always very small in the advective 
flow compared to that of a Keplerian disk near the outer boundary of 
the advective region. Thus some improvements of the disk model at the 
transition region is needed. Since major reactions are closer to the black 
hole, we believe that such modifications of the model would not change 
our conclusions. The neutrino luminosity in a steady disk is generally very 
small compared to the photon luminosity [6], but occasionally, it is seen to 
be very high. In these cases, we predict that the disk would be unstable. 
Neutrino luminosity from a cool advective disk is low. 

In all the cases, even when the nuclear composition changes are not 
very significant, we note that the nuclear energy release due to exothermic 
reactions or absorption of energy due to endothermic reactions is of the 
same order as the gravitational binding energy release. Like the energy 
release due to viscous processes, nuclear energy release strongly depends 
on temperatures. This additional energy source or sink may destabilize the 
flow [6]. 
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